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Abstract 

Based on the relationship between thermodynamics and gravity we propose, with the aid of 
Verlinde’s formalism, an alternative interpretation of the dynamical evolution of the Friedmann- 
Robertson-Walker Universe. This description takes into account the entropy and temperature 
intrinsic to the horizon of the universe due to the information holographically stored there through 
non-gaussian statistical theories proposed by Tsallis and Kaniadakis. The effect of these non- 
gaussian statistics in the cosmological context is change the strength of the gravitational constant. 
In this paper, we consider the tcCDM model modified by the non-gaussian statistics and investigate 
the compatibility of these non-gaussian modification with the cosmological observations. In order 
to analyze in which extend the cosmological data constrain these non-extensive statistics, we use 
type la supernovae, baryon acoustic oscillations, Hubble expansion rate function and the linear 
growth of matter density perturbations data. 
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I. INTRODUCTION 


There are theoretical evidences that the understanding of gravity has been greatly bene¬ 
fited from a possible connection to thermodynamics. Pioneering works of Bekenstein [l| and 


Hawking [2[ have described this issue. For example, quantities as area and mass of black 
holes are associated with entropy and temperature respectively. Working on this subject, 
Jacobson interpreted Einstein held equations as a thermodynamic identity. Padmanab- 
han ^ gave an interpretation of gravity as an equipartition theorem. Recently, Verlinde 
brought an heuristic derivation of gravity, both Newtonian and relativistic, at least for static 
spacetime. The equipartition law of energy has also played an important role. The analysis 
of the dynamics of an inhationary universe ruled out by the entropic gravity concept was 
investigated in . On the other hand, one can ask: what is the point of view of gravitational 
models coupled with thermostatistical theories and vice-versa? 


The concept introduced by Verlinde is analogous to Jacobson’s [3| one, who proposed a 
thermodynamic derivation of Einstein’s equations. The result has shown that the gravitation 
law derived by Newton can be interpreted as an entropic force originated by perturbations 
in the information “manifold” caused by the motion of a massive body when it moves away 
from the holographic screen. An holographic screen can be understood as a storage device 
for information which is constituted by bits. Bits are the smallest units of information. 
Verlinde used this idea together with the Unruh result and he obtained Newton’s second 
law. The idea of a entropic gravity/cosmology has been extensively investigated in different 
contexts, see [^, for recent results 

Moreover, assuming the holographic principle together with the equipartition law of en¬ 
ergy, the Newton law of gravitation could be derived. The connection between nonextensive 
statistical theory and the entropic gravity models 1^, 1I| make us to realize an arguably 
bridge between nonextensivity and gravity theories. In this paper we propose an alterna¬ 
tive interpretation of the dynamical evolution of the Friedmann-Robertson-Walker Universe 
(FRW) through non-gaussian statistical theories. We use the most recent observational data 
of Supernovae of Type la (SN la), Baryon Acoustic Oscillation (BAO), Hubble parameter 
and the growth function to investigate the cosmological consequences of such modifications 
through dark energy (DE) models. 

This paper is organized as follows. In section H we will make a brief review of the 
formulations concerning the non-gaussian of Tsallis and Kaniadakis statistics. In section 
HI we will present the formalism of Verlinde and its consequences for the gravitational 
framework. Section VI we introduce the modified dynamic FRW universe. In section V we 
use SN la, BAO, H{z) and f{z) data to constrain the non-gaussian statistics modifications on 
the tcCDM modified model. Lastly, Sec. VI briefly delivers our main conclusions and offers 
some hnal remarks. As usual, a zero subscript means the present value of the corresponding 
quantity. 


2 









II. NON-GAUSSIAN STATISTICS 


The objective of this section is to provide the reader with the main tools that will be 
nsed in the next sections. Althongh both formalisms are well known in the literatnre, these 
brief reviews can emphasize precisely that there is a connection between both ideas and that 
it was established recently 1^. The study of entropy has been an interesting task through 
recent years thanks to the fact that it can be understood as a measure of information loss 
concerning the microscopic degrees of freedom of a physical system, when describing it in 
terms of macroscopic variables. Appearing in different scenarios, we can conclude that 
entropy can be considered as a consequence of the gravitational framework [l|, [^. These 
issues motivated some of us to consider other alternatives to the standard Boltzmann-Gibbs 
(BG) theory in order to work with Verlinde’s ideas together with other subjects [ll| . 


A. Tsallis statistics 


An important formulation of a nonextensive (NE) BG thermostatistics has been proposed 
by Tsallis 12| in which the entropy is given by the formulation 


o'? 


Sq = ks- 


w 


q-1 


1 


i=l 


( 1 ) 


where pi is the probability of the system to be in a microstate, W is the total number of 
conhgurations and q, known in the current literature as Tsallis parameter or the nonextensive 
parameter, is a real parameter quantifying the degree of nonextensivity. The dehnition of 
entropy ([T]) has, as motivation, to analyze multifractals systems and it also possesses the 
usual properties of positivity, equiprobability, concavity and irreversibility. It is important 
to note that Tsallis’ formalism contains the BG statistics as a particular case in the limit 
g —)■ 1 where the usual additivity of entropy is recovered. Plastino and Lima [l3| used a 
generalized velocity distribution for free particles 141 


fq{v) =Bq 1 - (1 - g) 


mv 


1 

l-q 


( 2 ) 


2kBT\ ’ 

where Bq is a dependent normalization constant, m and v is the mass and velocity of the 
particle, respectively. They have derived a nonextensive equipartition law of energy whose 
expression is given by 


E. = 


-NknT 


5-3g‘ 

where the range of g is 0 < g < 5/3. For g = 5/3 (critical value) the expression of the 
equipartition law of energy, Eq. ([3]), diverges. It is easy to observe that for g = 1, the 
classical equipartition theorem for each microscopic degrees of freedom is recovered. 
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B. Kaniadakis statistics 


Kaniadakis statistics 15[, also called ^-statistics, similarly to the TT formalism general¬ 
izes the standard BG statistics initially by the introduction of ^-exponential and K-logarithm 
dehned by 


ln«(/) = 

where the following operation being satished 


exp,,(/) = (y/f + , 

r - n 


2k 


ln^(exp^(/)) = exp^(ln^(/)) = / . 


(4) 

(5) 

( 6 ) 


By Eqs. (jl]) and ([5]) we can observe that the ^-parameter deforms the usual dehnitions 
of the exponential and logarithm functions. 

The K-entropy associated with this K-framework is given by 


S.if) = - / d^Pf 




2k 


(7) 


which recovers the BG entropy in the limit k —>■ 0. It is important to mention here that the 
K-entropy satished the properties of concavity, additivity and extensivity. Tsallis’ entropy 
satishes the property of concavity and extensivity but not additivity. This property is not 
fundamental, in principle. The ^-statistics has been successfully applied in many exp erimen- 
tal fronts. As examples we can mention cosmic rays 1^, quark-gluon plasma 17(|, kinetic 
models describing a gas of interacting atoms and photons 1^ and hnancial models 19|. 

The kinetic foundations of ^-statistics lead to a velocity distribution for free particles 
given by 201 


Mv)=\ + 


'2kBT' 


2kBT 


( 8 ) 


The expectation value of is given by 


/ „,2 > _ f-^^dv 




(9) 


Using the integral relation 21| 


dxx'' ^exp^(—x) = 


2 

K 

r r(ifo-i) 

1 -|- r 



r(r), 


( 10 ) 


we have that 
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( 11 ) 


< v‘^ >« = 


2/cbT 1 1 + 1/t r(^ -1) r(^ +1) 
m 2Kl + lKT{± + l)T{±-iy 


The K-equipartition theorem is then obtained as 


E = V-i±i^hinl)£<i±i 

- 2 2Ki + |Kr(i + |)r(i ■ 


keT. 


( 12 ) 

2"' y 2 k , ' 4^ y 2 k 4-^ 

The range of^isO < k < 2/3. For k = 2/3 (critical valne) the expression of the 
eqnipartition law of energy, Eq. (Ha, diverges. For k = 0, the classical eqnipartition 
theorem for each microscopic degrees of freedom can be recovered. 


III. VERLINDE’S FORMALISM AND THE MODIFIED GRAVITATIONAL CON¬ 
STANT 


The formalism proposed by E. Verlinde derives the gravitational acceleration by nsing, 
basically, the holographic principle and the eqnipartition law of energy. This model considers 
a spherical snrface as the holographic screen, with a particle of mass M positioned in its 
center. A holographic screen can be imagined as a storage device for information. The 
nnmber of bits (the term bit means the smallest nnit of information in the holographic 
screen) is assnmed to be proportional to the area A of the holographic screen 

(13) 

T 

where A = and Ip = ■>/ Gh/c^. In Verlinde’s formalism we assnme that the total energy 
of the bits on the screen is given by the eqnipartition law of energy 

E = ]^NkBT. (14) 

It is important to mention here that the usual eqnipartition theorem, Eq. (ITT)) , is de¬ 
rived from the usual BG thermostatistics. In a nonextensive thermostatistics scenario, the 
eqnipartition law of energy will be modified in a sense that a nonextensive parameter q 
will be introduced in its expression. Considering that the energy of the particle inside the 
holographic screen is equally divided through all bits then we can write the equation 


Mc^ = -NkeT. 
2 


Using Eq. fll3p and the Unruh temperature formula 

, ^ 1 ha 

ksT = - -, 

27r c 

we are in a position to derive the (absolute) gravitational acceleration formula 


(15) 


( 16 ) 
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_ _ M 

h 


( 17 ) 


We can observe that from Eq. fll7l) the Newton constant G is just written in terms of the 
fundamental constants, G = IpC^/h. 

As an application of NE equipartition theorem in Verlinde’s formalism we can use the NE 
equipartition formula, i.e., Eq. ([3]). Hence, we can obtain a modified acceleration formula 
given by 11 | 


where Gq is an effective gravitational constant which is written as 

G. = l^G. 


(18) 


(19) 


From result (IT^ we can observe that the effective gravitational constant depends on the 
NE parameter q. For example, when g = 1 we have Gq = G (BG scenario) and for q = 5/3 we 
have the curious and hypothetical result which is Gg = 0. This result shows us that q = 5/3 
is an upper bound limit when we are dealing with the holographic screen. Notice that this 
approach is different from the one demonstrated in 2^, where the authors considered in 
their model that the number of states is proportional to the volume and not to the area of 
the holographic screen. 

On the other hand, if we use the Kaniadakis equipartition theorem, Eq. (IT^ . in Verlinde’s 
formalism, the modihed acceleration formula is given by 


G« = G, 


M 


where G^ is an effective gravitational constant which is written as 


Gk — 2k- 


i + |Kr(i + |)r(i 


i) 


G. 


( 20 ) 


( 21 ) 


From result fl2T|) we can observe that the effective gravitational constant depends on the 
K parameter. For example, from the Gamma functions properties we have that 


lim 


r(^- 


2 ) r(— +1) 

\2k ^ 


= 2k. 


^0 rf— + -j rf— - -) 

\2k ^ \2k 4 / 

Then, using Eq. fl2ip we obtain for k = 0 that G^ = G (BG scenario). 


( 22 ) 
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IV. DARK ENERGY MODELS THROUGH NON-GAUSSIAN STATISTICS 


It was demonstrated in ll(] that one modification in the dynamics of the FRW universe 
in NE Tsallis’ statistics can be obtained simply by making the prescription G ^ Gq = 
(5 — 3q)G/2 in the standard field equations. From this proposal in new cosmological 


constraints on the parameter q were obtained. Analogouly to Tsallis’ statistics, we can 
modify the Friedmann’s equation in the Kaniadakis framework by making the prescription 
G —)■ G^, with Gk given by fl^ . Thus, for a homogeneous and isotropic universe filled by 
perfect fluids, the equations of motion for non-gaussian statistics can be written as 


and 


tt 2 k 



(23) 


H + — —^Gq(K)(p + 3p), (24) 

where FT = a/a is the Hubble function, Gq(K) denotes the effective gravitational constant 
in the Tsallis (Kaniadakis) formalism and p and p are, respectively, the total density and 
pressure of the fluid. These equations can be combined to obtain the conservation equation. 


p + 3iF(l + '«;)p = 0. (25) 

For a FRW universe pervaded by radiation, non-relativistic matter (baryonic plus dark 
matter) and some sort of dark energy to take into account the late time cosmic acceleration, 
the Friedmann equation fl2^ becomes 

^2 ^ ~ ^ ^ ^ + ^xof (o)^ (26) 

where = ^'^Gpio /is the density parameter of the i-th component {i = 7 , m, and a; 
for radiation, matter (baryonic more dark) and dark energy, respectively), Qk,o = ~k/HQ is 
the curvature density parameter and 


\ Px -3 

/ (a) = — = a exp 
PxO 


-3 


w{a')da' 


(27) 


is the density ratio for a dark energy fluid with a generic equation of state parameter w{a) = 
Px/Px- In the above equations the subscript 0 denotes the observable at present time. 

In the non-gaussian scenario, the normalization condition reads as 


^ 7,0 + f^m ,0 + ^fc ,0 + — 77 ;-) (28) 

which is an interesting result since it can show us, one more time, that the value q = 1 
{k = 0) recovers the standard normalization condition. Values of q > 5/3 and k < 0 which 
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brings a negative normalization condition makes no sense. 


In this paper we assume spatial flatness and we restrict ourselves to the case w = const 
so that Eq. fl26|l becomes 




G 


q{K) 


HI 


G 


qO ^ hi. 


^ + 


G 


G, 


q(K) 


Q H 


mO 


(1 + 


w) 


(29) 


From Eqs. flT^ . fl2l]) and fl28|l it is possible to note that the parameters q and n affects the 
energy balance of the universe. If q is greater (smaller) than one, the effective gravitational 
constant is smaller (grater) than G so that more (less) dark energy will be required to provide 
the observed late time universe acceleration. By its turn, the Kaniadakis framework, if k > 0, 
the gravitational field is weaker than in the gravitational held in the standard BG scenario 
so that we will need more dark energy to accommodate the cosmic acceleration. Since k > 0, 
the Kaniadakis statistics is more restrictive than Tsallis statistics. In the next section, we use 
some of available cosmological observations to obtain new constraints on the non-gaussian 
statistical parameters q and k. 


V. OBSERVATIONAL CONSTRAINTS 


In order to constrain the parameters (q,K and w) we perform ajoint analysis involving 
the 580 SNe la distance measurements of the Union2.1 data set 2^, the 30 measurements 
of the Hubble parameter H{z) given in Table 4 of Ref. 2^, 17 measurements of the growth 
function f(z) listed in Table I and the six estimates of the BAG parameter given in Table 3 
of Ref. Ei. 


A. Type la Supernovae 

The SNe la sample of the Union 2.1 is given in terms of distance modulus /i. Theoretically, 
the distance modulus is given by 

fif"{zi) = 5\og Hodiizi) + 51og(3/h) +40, (30) 

where h = i7o/100 Km • s“^ • Mpc“^ and 

is the luminosity distance for a spatially flat universe. The usual function is calculated as 










( 32 ) 


2 ^ (fif-fif)" 

= 2^ -- 

i Mi 

where /i°^^ is the observed value of the distance modulus at redshift Zi and cr^, its uncertainty. 
In our analysis we treat the Hubble constant Hq as a nuisance parameter and marginalize 
over Hq so that the function to be minimized is 

XsN = ^ ^) (33) 

where the quantities A, B, C are given by: 

21 = (f - l‘t‘) (C'skh (4 - 4‘) . (34) 

B = Y.^c-4h (4 - 4‘) (35) 

i 

and 

c = E(Cst).i (36) 

n 

and {Cglf)ij is the inverse convariance matrix. 


B. Growth function 

In the linear regime the matter density perturbations 5 = 5pm/Pm satisfies 

h + 2H6 - AirGeffPruS = 0, (37) 

where Geff is the effective gravitational constant for a given theory of gravity. For the models 
studied in this paper, G^jf is given by flT^ for Tsallis statistics and by (1^ for Kaniadakis 
statistics. By dehning the growth factor / = d In h/ In a, this second order time differential 
equation is reduced to 

/' + y+(2-|n.„)/-|n^" = o, (38) 

where f = df/dlna and 


ntif = 


STrGeffPv 

37/2 




-3 


^ 7 , 0 ® ^ A AtnnflCl ^ -|- 


G 


G, 


q (k) 


^7,0 ^ 


mO 


q— 3 { 1 +w) 


( 39 ) 


In Table I we list the 20 measurements of /. The usual function is calculated as 
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Table I: Currently available data for growth rates used here. 


z 

/ 


Ref. 

0.02 

0.482 

0.09 

[27] 

0.067 

0.56 

0.11 

[28] 

0.11 

0.54 

0.21 

[29] 

0.15 

0.49 

0.14 

[30] 

0.15 

0.51 

0.11 

[31J; [32] 

0.22 

0.60 

0.10 

[33] 

0.32 

0.654 

0.18 

[3^ 

0.34 

0.64 

0.09 

[3^ 

0.35 

0.70 

0.18 

m 

0.41 

0.70 

007 

[33] 

0.42 

0.73 

0.09 

\37] 

0.55 

0.75 

0.18 

[38] 

0.59 

0.75 

0.09 

[32 

0.60 

0.73 

0.07 

m 

0.77 

0.91 

0.36 

[30] 

0.78 

0.70 

0.08 

[33] 

1.4 

0.90 

0.24 

[39] 

2.125 

0.78 

0.24 

m 

2.72 

0.78 

0.24 

m 

3.0 

0.99 

0.24 

[29] 




20 


_ j^theo'^2 


(Jf 


(40) 

i=l fi 

where is the observed value of the growth function at redshift Zi, af. its uncertainty and 

jtheo value of f{zi) provided theoretically. In order to obtain we use the approx¬ 
imation f{z) ~ 2 ^, where 7 is the growth index and depends of the underlying 

cosmological model. For the tcCDM model, 7 = 3(w — l)/( 6 tc — 5). 


C. Hubble parameter 


In Ref. j4l|, Jimenez has developed a method to use the relative age of old and passive 
galaxies, dz/dt, to infer the Hubble parameter as a function of the redshift. 


H{z) = - = 
a 


1 dz 


1 + z dt 


(41) 
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The cosmic chronometers method to measure H[z) does not depend of any integrated dis¬ 
tance measurement over redshifts and is independent of cosmological models. We use 30 
H{z) data obtained from the cosmic chronometers method listed in Table 4 of 2^. The 
usual function is calculated as 


30 


/ zjobs u 

xi = E-' ' 


theo\2 


i=l 


a 


(42) 


Hi 


where is the observed value of the Hubble parameter at redshift Zi, a Hi its uncertainty 
and the value of the Hubble parameter at Zi provided theoretically. 

D. BAO 


Our BAO analysis is based on the BAO parameter: 


A(z) = Dv(z)^n'S^Hl 


where 


Dy(z) = 


dz' \ 2-1 1/3 


H{z)\zJa H(z') 


(43) 


(44) 


is the so-called dilation scale and = Gq{k)Vtrafl is the matter density modihed to take 
into account the effects of the non-gaussian statistics. Here we use the six estimates of the 
BAO parameter given in Table 3 of Ref. 2^ . The usual function is calculated as 


Xbao ~ 


E 


^j^obs _ j^theo'^z 


Hi 


(45) 


4=1 

where is the observed value of the BAO parameter at redshift Zi, CAi its uncertainty 
and the theoretical value of the BAO parameter at redshift Zi. 

As the likelihood function is dehned by L oc exp(—x^/2) the said values follow from 
minimize the quantity Xtotai = + x]r + X/ + Xbao- following we use = 0.27. 


E. Results 


Figured] shows the 68% and 95% confidence regions in the q — w plan (left panel) and in 
the K — w plan (right panel) for the modified wCDM model. At 2ct, the best fit points are: 
(q = 0.981°;°^^, w = -1.02l°;°si) with xLn = 588.32 and (k = 0.00^°;^^°, w = -1.01;°;°^^) 
with Xmin = 589.25. Quintessential {w > —1) and phantom {w < —1) fluids are compatible 
with both, Tsallis and Kaniadakis statistics. As we can see, the cosmological observations 
used in this work are completely compatible with non-gaussian statistics. However, it is 
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1.1 




Figure 1: la and 2cj confidence regions for the modified rcCDM model from a joint analysis with 
SNIa + f + ff(z) +BAO data. Left: Tsallis’ statistics. Right: Kaniadakis Statistics. 

important to note that these data does not exclude the BG statistics, i. e., g = 1 and k = 0. 
Tsallis’ statistics provide a better fit to the data than Kaniadakis. For Tsallis’ statistics 
0 < Geff/G < 2.5 while for Kaniadakis’ statistics 0 < Geff/G < 1. Thus, we can to 
attribute this result to the additional freedom allowed by Tsallis’ statistics since the data 
seems to choose a slightly greater gravitational constant. 


VI. CONCLUSION 

Currently, two extensions of standard statistical mechanics, known as Kaniadakis and 
Tsallis statistics, has been used to explain a very large class of phenomena observed experi¬ 
mentally in different areas, e.g, in low and high energy physics, astrophysics, econophysics, 
biology, ect. In this paper, we have explored the possibility of one modihcation in the dy¬ 
namics of the FRW universe obtained through an entropic force theory generalized for the 
Kaniadakis and Tsallis statistics. 

From the combination SNIa-|-BAO-|-Ff(;2) + f{z) datasets we have obtained new cosmo¬ 
logical constraints over q and k parameters. 

To sum up, in this paper, we have obtained new values for the free parameters that 
characterizes the non-gaussian statistical theory proposed by Kaniadakis and Tsallis.Based 
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in the data used in this paper, we note that non-gaussian statistics can not be rule out by 
cosmological observations although the BG statistics remains in fully agreement with the 
data. From the results obtained here we can conclude that the parameters q and k affects 
the energy balance of the universe. The gravitational held is more weak for q > 1 and 
K > 0 so that we need more dark energy than we would have if we consider the standard 
BG scenario. For 0 << 1 the strength of the gravitational held greater than in the standard 
BG scenario and less dark energy is required to explain the cosmological observations. The 
results obtained in this paper favor a slightly strong gravitational held. 
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